A decomposition of light’s spin angular momentum density

Light carries intrinsic spin angular momentum (SAM) when the electric or magnetic field vector rotates over time. A familiar vector equation calculates the direction of light’s SAM density using the right-hand rule with reference to the electric and magnetic polarisation ellipses. Using Maxwell’s equations, this vector equation can be decomposed into a sum of two distinct terms, akin to the well-known Poynting vector decomposition into orbital and spin currents. We present the first general study of this spin decomposition, showing that the two terms, which we call canonical and Poynting spin, are chiral analogies to the canonical and spin momenta of light in its interaction with matter. Like canonical momentum, canonical spin is directly measurable. Both canonical and Poynting spin incorporate spatial variation of the electric and magnetic fields and are influenced by optical vortices. The decomposition allows us to show that a linearly polarised vortex beam, which has no total SAM, can nevertheless exert longitudinal chiral pressure due to equal and opposite canonical and Poynting spins.


Introduction
Optical angular momentum is contributed to by physical orbital and spin parts [1][2][3][4].The orbital component can be either extrinsic or intrinsic and is associated with the spatial structure of the light field.Orbital angular momentum (OAM) is most famously stirred into light by optical vortices, which have helical phase fronts twisted around a onedimensional line singularity [5][6][7][8].Spin angular momentum (SAM), meanwhile, develops as the electric and magnetic field vectors rotate during oscillation [9], and is an intrinsic quantity, independent of co-ordinate origin.Both OAM and SAM arise from chiral structures in light, and as such, can access the chirality of matter, enhancing, pushing, twisting, and torquing with or without a parity bias [10,11].Many commonly observed chiroptical interactions, such as circular dichroism, arise from the coupling of chiral matter to a photon's spin state.Recent attention has also turned to OAM-dependent photon absorption and scattering, with a range of effects theorised [12][13][14][15].Chiroptical effects involving OAM often hinge on the strong longitudinal field component introduced by beam focussing [16,17].Unlike for photons with opposite spins, however, preferential absorption of photons with different OAM handedness does not occur in the dipole approximation for paraxial light [18,19].
A well-known orbital and spin separation can be performed locally in the Poynting vector [20] (kinetic momentum density, when divided by c 2 ).The Poynting vector may be expressed as a sum of the orbital current, pointing in the direction of canonical momentum, and spin current, proportional to the virtual spin momentum.This is a physically meaningful decomposition, tying in to the famous Abraham-Minkowski dilemma [21,22].
Surprisingly to most readers, a similar vector decomposition can be performed on the electromagnetic spin density S, which is split into, perhaps confusingly, orbital-like and spin-like contributions to the total electromagnetic spin.We name these two terms canonical and Poynting spin.This little-known decomposition is the focus of our work-amidst discussion by Shi et al who, without giving the complete decomposition explicitly, linked the two terms to longitudinal and transverse spin in some specific cases [23][24][25], no clear picture exists in the literature on the physical significance of the two terms in general, in the same way as for the decomposed Poynting vector.The decomposition as we present has been expressed previously (e.g.[26]), but with almost no discussion on the function of terms as component parts of S.
We therefore have two objectives: to clarify and understand the meaning of the orbit-like and spin-like components of S, and to emphasise their value for the research community in interpreting electromagnetic SAM density.The component parts of light's total spin angular momentum density, visualised at a single point (red circle) in the 3D interference of monochromatic plane waves.Enlarged electric and magnetic field ellipses and coloured vector arrows, representing the total spin S, electric and magnetic spins S e and S m , and the canonical s c and Poynting spins s p of the decomposition Eq. ( 4), are plotted.The electric and magnetic spin vectors are normal to the E and H polarisation ellipses, according to the right-hand rule and the sense of rotation of the instantaneous field vectors (see separated diagrams below the main combined image).Projected on the xy plane are the Poynting vector streamlines, which in this case organise a clear curling structure.This provides a visual aid for Poynting spin (large black arrow), which is informed by three-dimensional curling of the Poynting vector.In structured light where the E and H ellipses are de-coupled, the canonical spin and therefore chiral pressure (Eq.( 5)) can point in a completely different direction to the any of the total spin or its electric and magnetic contributions.
For example, by use of the Maxwell stress tensor, one understands that the canonical and Poynting spin vectors are responsible for different interactions between light and chiral particles.While the two decomposed spin vectors are organised independently in general monochromatic light, a remarkable implication lies in beams carrying OAM, such as vortex beams.OAM imparts a longitudinal component to both canonical and Poynting spins, even in the absence of longitudinal total spin.By applying the decomposition to a linearly polarised vortex beam, we show that a total spin-free longitudinal chiral force exists-a force which can act even under the dipole approximation, despite originating from the beam's OAM.
We generalise the spin decomposition to time-dependent fields and four-vector representation, unlocking a deeper understanding of the behaviour of the two terms.We also show that the ability to split spin into two terms appears to be a general feature of wave fields, similarly to the decomposition of the Poynting vector, including those of the linearised theories of gravity [27] and acoustics, which have distinctly different vector structures.Throughout, we use bold latin letters E and H to represent time-harmonic complex electromagnetic field phasors, which are a function of position only [e.g.E(r) = E 0 exp(ik • r)], and scripted characters E and H denoting real, time-varying field vectors (e.g.E(r, t) = Re{E(r) exp(−iωt)}).

Electromagnetic spin decomposition
We will first contextualise the decomposition of electromagnetic SAM density by recalling a well-known orbital and spin decomposition which exists in free-space and time-harmonic fields for the time-averaged Poynting vector, which represents the flux of active power in the light field.The Poynting vector is separated into the orbital and spin currents [20,22], defined using the complex phasors E and H.The inner product notation is a • (∇)b = a x ∇b x + a y ∇b y + a z ∇b z .When divided by c 2 , P has units of a momentum density and is termed the kinetic momentum density, while the orbital current p o becomes the canonical momentum density, corresponding to the expectation value of linear momentum carried by photons at each point in space (unlike kinetic momentum, the canonical momentum is a directly measurable quantity).The net flow of the instantaneous Poynting vector into or out of a volume relates to the change in electromagnetic energy density over time, according to a continuity equation.An adjacent conserved quantity of light is its chirality, which has its own continuity equation and an associated chirality flux [28][29][30].Chirality χ is one of an infinite hierarchy of conserved quantities in linear, non-dispersive media, to which helicity h also belongs [9,30]-helicity and chirality are often conflated because they (and their fluxes) are proportional by a factor of ω 2 in monochromatic waves, though in general they are distinct.In monochromatic light the flux of helicity corresponds to the time-averaged SAM density (often referred to simply as 'spin'), given by, Spin's individual electric and magnetic contributions point in the normal direction to the electric and magnetic polarisation ellipses, drawn over time by the instantaneous vectors E(r, t) and H(r, t).Like the Poynting vector, and using the same procedure, spin S may be split into a sum of two vectors: the canonical spin s c and the Poynting spin s p .This decomposition is our main focus and is, explicitly, Any subsequent mention of 'canonical spin' in this work always refers to the first term, s c = 1 4ω 2 Re{E * • (∇)H − H * • (∇)E}, and 'Poynting spin' to the second, s p = 1 2ω 2 ∇ × P. Their sum S (Eq.( 3)) we refer to as total spin.This naming scheme reflects that of the decomposed kinetic momentum, split into canonical momentum and spin momentum.

Physical interpretation
Previous works [23,24] have linked the canonical s c and Poynting spin s p vectors to the longitudinal and transverse spin of light, respectively.However, this is only true in a very limited number of cases such as linearly polarised evanescent waves.It cannot be a general result because well-defined longitudinal and transverse directions do not exist for many-plane wave or multiple beam interference.In fact, it does not hold even when there is a well-defined longitudinal direction: a circularly polarised evanescent wave has both longitudinal and transverse spin components contained in the Poynting spin vector.
One way to gain an understanding of the physical significance of the canonical and Poynting spins is to study the interactions between light and chiral matter.Chiral light can exert preferential forces on enantiomers with opposite handedness, with direct proportionality to the light's SAM.When light shines on a particle much smaller than the wavelength (Rayleigh regime) it gets polarised and acquires an electric p and a magnetic m dipole.In the linear regime, the dipole moments are proportional to the incident fields, p = α e εE + iα c H/c and m = α m H − iα c E/η, where α e , α m and α c are the electric, magnetic, and chiral polarisabilities of the particle.The chiral polarisability α c is a pseudoscalar so it changes sign between the two enantiomers (mirror-reflected versions), and vanishes unless the matter is chiral.Under these assumptions, the illuminating light exerts a chiral optical force on the particle (which changes sign for the different enantiomers) given by [31][32][33][34][35]: where k = ω/c is the wavenumber, h = −Im(E * • H)/(2ωc) is the cycle-averaged optical helicity density [17,36], and S e and S m are the electric and magnetic parts of the total spin S from Eq. (3).One immediately sees that the canonical spin s c , one of the two terms in the spin decomposition [Eq.( 4)], appears in this force equation and is directly responsible for the chiral pressure (in this context, the product ks c is sometimes called the chiral momentum [26]), while the total spin appears in the spin recoil term.Poynting spin contributes (together with canoncial spin) only to the higher-order recoil force caused by the unbalanced radiation pattern of an electric-magnetic dipole.This also implies that in an electromagnetic field whose SAM is pure Poynting spin (i.e., zero canonical spin), the particle will experience no chiral pressure, only helicity gradients and relatively weaker spin recoil terms of the force.An alternative route to the physical meaning of the canonical and Poynting spins is found purely in the two terms' mathematical expressions as shown in Eq. ( 4).The Poynting spin, being proportional to the curl of the Poynting vector, is easily interpreted as the vorticity in energy flow.The expression for the canonical spin is much harder to interpret initially, but begins to unravel if we decompose a general electromagnetic field into circularly polarised plane waves.Any arbitrary electromagnetic field, no matter how complicated, can be expressed in momentum space by an angular spectrum (an infinite sum of plane waves of different wavevectors, weighted by an amplitude function).Using this property, we can further probe the canonical spin term in monochromatic light by separating electric and magnetic fields into two component fields of opposite helicities, indicated by the + and − subscripts, where ê± (k) represent the circularly polarised unit vectors for each plane wave with wave-vector k (see, e.g., [37]).
The magnetic field's associated helicity components are obtained from Faraday's law ∇ × E = iωµH as, where we used the property (k/k) × ê± = ∓iê ± .Equations ( 6) and (7) show that the helicity-separated electric and magnetic fields are related by H ± (r) = ∓iE ± (r)/η, characteristic of a circularly polarised plane wave-this is true not only in the spectral representation, but also in the spatial representation, for any arbitrary field (a property explored in depth in [38]).This allows us to substitute helicity-separated fields into many dual quantities, including canonical momentum and canonical spin, gaining further insight.Simply substituting E = E + + E − and H = −i (E + − E − ) /η into the expression for orbital current [p o of Eq. ( 2)], we find, after some algebra, showing that the helicity segregation of the E and H fields translates to a separation of contributions to the field's momentum p = p + + p − by photons of positive and negative helicity.Note that orbital current is proportional to canonical momentum by p o = c 2 p; here, p + and p − are helicity-separated momentum densities.Making the same substitution in the expression for canonical spin s c [first term of Eq. ( 4)] illuminates its meaning, Canonical spin is proportional to the difference in linear momentum densities carried by photons of oppositely signed helicity.This brings a clear physical interpretation of the canonical spin, and also helps explain why the chiral pressure force acts in the direction of s c , as photons of opposite helicities are absorbed or scattered in different amounts by chiral particles.In a general structured field, p + − p − ̸ = p + + p − and hence canonical spin is decoupled from any local longitudinal direction defined by canonical momentum.

Mathematical derivation
In this section, we briefly lay out the mathematical steps taken to arrive at Eq. ( 4), before developing the expression into a more fundamental 4-vector description.For each of the canonical and Poynting spin terms, the 4-vector decomposition incorporates a time component which characterises how the two terms transform differently between reference frames.

3-vector decomposition
Maxwell's equations enable both decompositions of the Poynting vector Eq. ( 2) and spin Eq. ( 4) by bringing the relevant vector from its usual representation into a form which can be separated into two terms using a vector identity, where a and b are arbitrary vectors, and the previously unseen notation on the right hand side means (a In time-harmonic fields, the phasors E and H can be substituted for curls of their counterpart via Faraday's and Ampere's laws.Substituting for the un-conjugated phasors in the spin vector definition Eq. ( 3) gives an expression in the form of Eq. ( 10), Applying the vector identity Eq. ( 10), we have, Gauss' law in free space (∇ • E = 0 and ∇ • H = 0), combined with a second vector identity , extracts the curl of the Poynting vector from the second Re{} term of Eq. ( 12), completing the final step in obtaining the electromagnetic spin decomposition as we present in Eq. ( 4).
In addition to the monochromatic case, canonical and Poynting spin analogies for the flow of chirality may also be defined in polychromatic fields, expressed using the full time-dependent fields E(r, t) and H(r, t) and stemming from the instantaneous flow of chirality F (r, t) [28][29][30]39].
Chiral flow F , unlike the flux of helicity in time-dependent fields, is defined without the vector potentials A and C, which avoids gauge complications and contains a direct relation to the instantaneous Poynting vector P = E × H.For monochromatic light, F is not time-dependent and becomes proportional to the field's SAM density S. Equation ( 14) is immediately presented in a form which can be broken into two vectors using the two identities detailed in this section, giving the extension of the time-harmonic spin decomposition to time-dependent fields, The above is a more general decomposition of chirality flow than Eq. ( 4), and is valid for polychromatic aperiodic light, at every instant in time.It is worth stressing again that the flow of chirality F is a different quantity to the flow of helicity (SAM density) for general time-dependent fields.To obtain the helicity flow equivalent of Eq. ( 15) in the Coulomb gauge, one simply substitutes E → A and H → C. Equation ( 15) assumed a source-free medium, such that ∇ • E = ρ/ϵ = 0. Dropping this assumption, Eq. ( 15) can be further generalised by adding a third term, ρ ϵ H , to the decomposition.

4-vector representation
In an attempt to extract new physical insight, we reproduce the above spin decomposition Eq. (4) in 4-vector notation.This notation is particularly elegant and useful to verify the Lorentz covariance of quantities and equations.
We will begin by introducing some common definitions.Electrodynamics is intrinsically relativistic, and several related physical quantities can be expressed via well-known 4-vectors, such as the 4-potential (grouping the scalar and vector potentials) and the 4-current (grouping the charge and current densities).4-vector notation is well known to provide an efficient way to formulate Maxwell's equations.Indeed, the time-dependent Maxwell's equations (using scripted vectors) in vacuum are a combination of four equations, two scalar equations ∇ • E = 0, ∇ • H = 0 and two vector equations In 4-vector notation, however, we can write Maxwell's equations as: Throughout, the greek indices µ, ν, α, β... run from 0 to 3, where the 0 th component labels the time direction.Repeated indices are summed over.Note that the distinction between subscript/superscript placement of indices are important, as A µ = η µν A ν where the so-called Minkowski metric η µν = diag(−1, 1, 1, 1).We also use Roman indices i = 1, 2, 3 to denote the three-dimensional space.In Eq. ( 16), ϵ αβγλ is the Levi-Civita symbol [40], while F αβ is the field strength tensor that conveniently packages the electric and magnetic fields: The field strength tensor can be expressed as , where A α is the 4-vector potential, that encases the scalar potential φ and the vector potential A and ∂ µ is a generalisation of gradient: Equation ( 16) is a striking example of the neatness of the tensorial notation: the physical content is the same, but the tensorial notation is clearer and exhibits Lorentz covariance.For time-harmonic fields, these tensors will also have a phasor representation (e.g.F µν (r, t) = Re{F µν (r) exp(−iωt)}), including the time components of the gradient, which for phasors becomes ∂ 0 = −iω/c.We now have the tools to write, using phasors, the time-averaged spin decomposition Eq. ( 4) in tensorial notation.The first step is to equivalently express the spin, defined in Eq. ( 3), using spatial indices to replicate the vector operations, Then, the equivalent of the proposed decomposition in Eq. ( 4) is: The next step is to find a 4-vector, for which the spatial part would reduce to Eq. (20).It turns out that this 4-vector quantity is exactly the helicity density and flux, a 4-current density associated with conserved helicity [41] where C µ = (ψ, −iωµ −1 0 (∇ × A)) ⊺ is the magnetic equivalent to the 4-potential A µ , and G µν = ∂ ν C µ − ∂ µ C ν is the corresponding field strength tensor.Note that Eq. ( 21) is time-averaged and the fields are in their phasor representation.We include the instantaneous variant in the supplementary information.Working in the Coulomb gauge, i.e., the scalar and vector potentials are chosen such that ∇ • A = ϕ = 0 and ∇ • C = ψ = 0, and after some algebra (see supplementary Mathematica file [42]), we find that the spatial part of the 4-vector S µ reproduces Eq. ( 20), whilst the temporal part is the cycle-averaged helicity density h = −Im{E * • H}/(2ωc), This is consistent with the fact that h and S are density and flux associated with integrated helicity, the conserved quantity associated with the dual symmetry [9,30] (hence ∇ • S = 0 for free-space monochromatic fields).This four-vector current S µ has a decomposition equivalent to Eq. ( 4), where the individual components are given by, reproducing the spin decomposition in 4-vector notation.Note that the temporal part of the total spin, the helicity density, is carried only by the canonical 4-spin term while the time component of the Poynting 4-spin is zero.
Let us briefly comment on the choice of the Coulomb gauge.As can be seen from the supplementary information, the decomposition of S µ is independent of the choice of frame and gauge.As already discussed in previous work (e.g.[9,41]), while the local density S µ is not gauge independent, upon integration one gets integrated helicity, which is a gauge-independent quantity that is conserved in any scenario where there is electromagnetic dual symmetry.Furthermore, the integral depends only on the transverse components of the potentials, which makes the Coulomb gauge the most convenient gauge to work with, and in this gauge, the helicity flux (spatial component of S µ ) also coincides with the optical spin density.Note, however, that the gauge fixed version of this decomposition is not invariant under Lorentz boosts, unless we simultaneously change our potentials to the Coulomb gauge associated with the new boosted frame.

Spin decomposition examples
It is instructive to apply the spin decomposition in Eq. ( 4) to some simple examples of light: a general evanescent wave and three commonly known focused beams, namely a Gaussian beam, a radial/azimuthal beam (l = 0), and an l = 1 vortex beam (all linearly polarised in the transverse plane).What differentiates each of these fields are their energy density and phase structures, due to one-dimensionally monotonic, doughnut or Gaussian real-space amplitude profiles, and the presence of OAM, two characteristics which reorient the Poynting vector throughout space and generate different amounts of Poynting spin.Poynting spin can emerge even when the instantaneous vectors E and H do not actually rotate, because its counterpart (canonical spin) is able to counterbalance the total SAM density of the field as necessary.Interplay between spin's canonical and Poynting components creates counter-intuitive effects, such as spin-free chiral interaction forces.

Evanescent wave
The electric and magnetic field phasors of an evanescent wave of arbitrary polarisation, propagating in the z direction (k z > k) and decaying along the x axis with decay constant γ = k 2 z − k 2 , are, where η = µ0 ϵ0 .Choosing complex values for A s and A p , which are the amplitudes of the evanescent wave's TE and TM modes, controls the wave's polarisation-a circularly polarised wave, for instance, has A p = ±iA s .The energy density of the wave decays in the x direction and is given by, Using our formulae, we can calculate the energy-normalised total spin of the evanescent wave, as well as its canonical and Poynting spins, The parameter ) is the degree of circular polarisation in the sense of a plane wave (σ = ±1 for circular polarisation, σ = 0 for linear polarisation).It is well-known that evanescent waves carry transverse spin independently of polarisation [26].This property is accounted for by the ŷ component of S, which is unaffected by the relationship between A s and A p and, interestingly, is a product solely of Poynting spin as was discovered in [23].Meanwhile, the evanescent wave acquires a longitudinal spin component if σ ̸ = 0, a component which is contributed to by both decomposed spins, and in different amounts.The over-generous canonical spin develops a larger ẑ component than is physical for the total spin of the wave (s c • ẑ > S • ẑ).Compensating, the Poynting spin's ẑ component points backwards to ensure (s c + s p ) • ẑ = S • ẑ (this becomes clear after substituting γ 2 = k 2 z − k 2 ).Decomposing an evanescent wave's spin reveals a physical distinction between its transverse and longitudinal spin components.Since Poynting spin is responsible for the transverse component of S, transverse chiral forces felt by an enantiomer arise as it recoils from its own radiation, rather than from a direct field interaction.This is consistent with the fact that perpendicular to the wavevector, there is no phase advance to twist the rotating field vectors into helices-the evanescent wave only acquires helicity if it carries longitudinal spin (h = W σ/ω). Canonical spin, purely longitudinal on the other hand, couples directly to chiral matter to produce a significantly stronger (relatively) preferential force.

Beams
Let us first consider a tightly focussed, linearly polarised Gaussian beam, whose energy density and polarisation in the focal plane is plotted in the top row of Fig. 2(a).In the transverse xy plane, the projection of the electric and magnetic field vectors are linear (x and y polarised respectively), though due to the tight focusing both E and H have significant longitudinal z components, such that they are elliptically polarised in 3D and contribute a circulation of transverse spin.The Poynting vector magnitude is not constant across the face of the beam and therefore has a non-zero curl in the transverse plane (except in the local maximum at the beam centre).This, the total spin and the decomposed canonical and Poynting spin of the Gaussian beam are plotted in Fig. 2(a)'s column (descending).Both E and H fields have z components and discretely symmetric polarisation ellipses in 3D (after a rotation of one field by π/2), which appears to reduce the canonical spin of the beam to zero, s c = 0.All transverse spin developed by focussing of a linearly polarised Gaussian beam stems from Poynting spin, S = s p .A similar conclusion could be made for an evanescent wave, making it tempting to argue for a general association between transverse spin and Poynting spin.This notion is not supported by the following example, however.
An l = 0 azimuthally polarised (in the electric field) beam and its spin decomposition is plotted in Fig 2(b).Considering dual spin (contributions from both E and H), the decomposition properties of this beam are identical for a radially polarised beam, which has an azimuthal magnetic field.Despite tight focussing, the azimuthal electric field vector does not carry a longitudinal component, as is well-known, and is completely linearly polarised in 3D.This means that the beam's wholly transverse spin is supplied by the magnetic field alone.Combining to give the total spin of the beam, both canonical and Poynting spins are non-zero-this contrasts with the Gaussian beam, whose electric and magnetic fields are both elliptical in 3D.Neither beam considered so far has twisted wavefronts carrying OAM, though due to their azimuthal transverse spin components, both the Gaussian and radial beam possess a certain chiral OAM.A chiral particle could access the spatial structure of the beam's spin field S, feeling an azimuthal force (clockwise or anticlockwise depending on the enantiomer) which would cause the particle to orbit the beam centre with no need for helical phasefronts.
Finally, we treat the x-polarised l = 1 vortex beam [Fig.2(c)], whose OAM leads to the most surprising spin decomposition features of the three.Numerically simulating such a non-paraxial, focussed beam is a challenging task because there is no unanimous 3D definition of a vector vortex beam (only scalar and paraxial beams are well-defined).The beams treated in this section are all electric linearly polarised and generated numerically using an angular spectrum integration technique [43] which, by convention, forces the electric field to exactly match the paraxial description of the beam in the focal plane.All unenforced field components (i.e., all magnetic field components and the longitudinal electric component throughout space, and the transverse electric components outside of the focal plane) are subsequently calculated via Maxwell's equations.While the beams are, therefore, physically valid, some magnetic-biased polarisation features develop in the x-polarised vortex beam considered here.In particular, although the electric field remains perfectly linear, the magnetic field alone gains a slight ellipticity in the transverse plane, contributing to a small longitudinal total spin outside of the beam centre as has been shown to exist in focussed linearly polarised beams [16].Longitudinal spin remains zero in the singularity, however.Breaking continuous rotational symmetry via x polarisation, both the vortex beam's electric and magnetic fields possess a z component from focussing, in contrast to the azimuthal beam for which the electric field is perfectly linearly polarised in 3D.The beam has helical wavefronts so that the average phase gradient (local wavevector) has an azimuthal component, also inherited by the Poynting vector (shown in the second row of Fig. 2(c); there is a small transverse component circulating in the xy plane).The Poynting spin s p ∝ ∇ × P, therefore, acquires a longitudinal component (even at the centre of the beam), in outright defiance of the fact that both transverse electric and magnetic fields are nearly linearly polarised, and zero in the centre, as is visible in the lowest plot of Fig. 2(c).To neutralise the longitudinal Poynting spin and suppress the total longitudinal spin, the z component of canonical spin s c must be non-zero too, satisfying s cz = −s pz for S z = s cz + s pz = 0 in the centre of the beam, and S z = s cz + s pz ≈ 0 elsewhere (see the plot second from bottom).In the middle of the beam, this counteracting canonical spin points in the negative z direction as a consequence of the beam's vortex handedness (right, topological charge +1), and would switch sign in a left-handed vortex.Returning now to the chiral force equation Eq. ( 5), we can infer that the beam's non-zero canonical spin, proportional to chiral momentum, produces a longitudinal chiral force, even though near to the vortex centre the E and H fields have negligible ellipticity in the transverse plane.Perhaps counter-intuitively, the longitudinal chiral force is strongest in the centre of the beam where the electromagnetic field is virtually zero (due, of course, to the longitudinal curl of the Poynting vector being maximal there).This is another example of spatially chiral light coupling to chiral matter which, even under the dipole approximation and in a linearly polarised field, can feel a discriminatory force by absorbing or scattering photons carrying OAM of a certain handedness.Only recently has OAM-dependent chirality, able to couple even to chiral dipoles, been demonstrated in focussed beams, resulting from strong longitudinal field components [17,44].We point out, however, that the chiral pressure force which we theorise at the centre of a linearly polarised vortex appears only to require that the Poynting vector has an azimuthal component (i.e., no requirement for longitudinal field components), and should still be present in a paraxial beam.Figure 2: Spin decomposition of non-paraxial beams.The beams, each of waist 1.5λ and separated in columns, are (a) a linearly polarised Gaussian beam (E T ||x), (b) an azimuthally polarised (E|| φ), l = 0 doughnut beam, and (c) a linearly polarised (E T ||x) vortex beam with topological charge l = 1.The top row of (isometric) plots across each sub figure column shows the beam energy density in colour, as well as the Poynting vector and electric (blue) and magnetic (green) polarisation ellipses, which are elliptical due to a significant z field component.Subsequent rows are vector plots of each beam's Poynting vector P, total spin S Eq. ( 3), canonical spin s c and Poynting spin s p from Eq. ( 4), respectively.White arrows are projections of the corresponding vector into the xy plane, while the red arrows are projections of the vector onto longitudinal yz and xz cut planes.Within each beam, arrows in the three spin decomposition plots are drawn to a consistent scale.Each non-paraxial beam is generated using an angular spectrum integration method [43].Defining the 3D vector vortex beam (c) is a difficult problem and the method we used produces a small and physical longitudinal spin (third row of (c)), contributed by the magnetic field, which would not be present in a (non-physical) perfect paraxial beam.

Analogies in other wave fields
Kinetic, canonical and spin momenta analogous to the terms in Eq. ( 2) can be identified more generally in other wave fields [45][46][47].This, and the ability to split the SAM density of light into two distinct terms poses another curiosity: what does a spin decomposition of another wave field look like, particularly if, unlike the electromagnetic field, its quanta are not spin-1?Treated in this section are acoustic and gravitational waves, which both depart from light's spin-1 structure.Despite their increased complexity, linearised theories exist for each of these fields under certain conditions.In a perfect fluid, acoustic waves are linear and oscillate longitudinally (spin-0), while gravitational waves are tensorial in nature (spin-2) and combine linearly when their amplitudes are sufficiently low.
In formulating equivalent expressions for total SAM, canonical and Poynting spins in linearised gravity, we use the Maxwellian representation of gravity in the weak field limit [27], i.e., waves propagating over flat spacetime at a large distance from their source.Gravitational waves detected on earth [48,49] arrive within this limit, and share some of light's characteristics; both fields have two polarisation degrees of freedom (in gravity these are the h + and h × polarisations) and are massless.Some of light's strangest behaviour, such as in evanescent fields, have also been predicted with additional properties in gravitational waves [50].Note that we take the expressions for the helicity and SAM density that are dual-symmetric (i.e., S = S e + S m rather than S = 2S e ), the same as in the case of electromagnetism, despite the fact that in the case of gravity we have no experimental evidence favouring its physical relevance over the asymmetric version [51].However, if we take the point of view of [9] that only the integrated helicity is a physically meaningful quantity, the two will be equivalent.
Table 1 summarises differences and similarities between linearised acoustics, electromagnetism, and linearised gravity with a focus on the spin-related quantities.An acoustic wave field can be described by a scalar pressure field P and a vector velocity field v which, in linearised acoustic theory, share a Maxwell-like relation [46,47].The derivation of the decomposition for the acoustic field is given in the supplementary material.A gravitational wave can be described using a metric perturbation h ij , which can be thought of as components of a three-by-three symmetric matrix.If we consider only i-th row or column of this matrix as a vector potential, A i = h ij êj , then the Maxwellian representation of gravity can be written in vector notation; the derivation of the spin decomposition then becomes simply Eq. ( 12).The acoustic helicity density is zero, which is a feature of spin-0 fields, but interestingly the canonical spin, which we showed in Eq. ( 41) to be related to helicity, is also zero.One can see that electromagnetic and gravitational waves are spin-1 and spin-2 respectively by taking a circularly polarised wave.The helicity density of a circularly polarised wave will be one(two) times the energy density and canonical spin will be one(two) times the canonical momentum for an electromagnetic(gravitational) wave.

Linearised Acoustics
Electromagnetism Linearised gravity Energy density  [45,46,50].For electromagnetism, the potential is considered to be in the Coulomb gauge.For linearised gravity, h ij are spatial components of the metric perturbation in the transverse-traceless gauge, êi are basis vectors, and any repeated indices are summed over (Einstein's convention).Parameters ϵ 0 = 1/(c 2 µ 0 ) = c 2 /(32πG) were chosen such that the time-averaged energy density takes the same form as the expression for the electromagnetic field.A larger version of this table using the index notation introduced in Sec.4.2 is given in the supplementary material.

Conclusions
We have discussed a known but overlooked decomposition which exists for the spin angular momentum density of light, in a similar way to the Poynting vector, which can be split into orbital and spin currents.Spin is decomposed into two terms, which we call the canonical spin s c and Poynting spin s p .We have further expressed the decomposition for time-varying fields, as well as in four-vector notation.The canonical spin is proportional to chiral linear momentum which can be transferred to chiral enantiomers in a positive or negative direction, depending on the enantionmer handedness.The resulting force is a chiral analogy to the (achiral) force due to radiation pressure-both are relatively strong first order forces and can act on matter in the dipole approximation.
The mechanism for preferential photon absorption or scattering by chiral matter is normally associated with the photon's spin state (i.e., circular dichroism), and less so to the sign of its OAM (although these interactions have received recent attention).We emphasised, however, that light is capable of exerting a first-order chiral force which does not directly depend on the total SAM of light, rather, its canonical spin, which can take a non-zero value even in linearly polarised fields.This is due to the second term in the spin decomposition, the Poynting spin, which by definition (being the curl of the Poynting vector) depends strongly on light's OAM: optical vortices twisting energy flow around the beam axis naturally generate non-zero Poynting spin, regardless of polarisation.We showed that under linear polarisation, an optical vortex's Poynting spin must be compensated by canonical spin, to produce a longitudinal chiral force in the absence of longitudinal SAM.Interestingly, this OAM-dependent chiral force is strongest in the centre of the vortex, where the electromagnetic energy density is minimal, along with achiral forces such as radiation pressure and gradient force.It is possible to trap atoms and molecules (using a blue-detuned wavelength with respect to a strong resonance) in dark spots [52][53][54].
The ability to decompose spin appears to be a general property of wave fields, and as we demonstrated, is straightforward to perform in linearised acoustics and gravity, whose quanta are considered spin-0 and spin-2 respectively.This fact could be of significant interest to a broader community, beyond optics.

Supplementary Information
This section supports the main text with four main components.In section 9.1, we provide an alternative version of Table I of the manuscript, expressed completely using index notation, which highlights the similarities and differences between vector and tensor waves.Secondly, section 9.2 is a discussion of the equivalent SAM density decomposition in linerised acoustics, which informs the first column of Table I of the manuscript.In section 9.3, details of the units of the quantities handled in the manuscript's 4-vector decomposition are given.Finally in section 9.4, the electromagnetic spin decomposition is presented in a general gauge, supporting comments made in the in 4-vector section of the main text.

Alternative table (index notation)
An alternative version of Table 1 from the manuscript is given below, where all quantities are expressed using index notation.

Decomposition in Linearised acoustics
An acoustic wave field can be described by a scalar pressure field P and a vector velocity field v which, in linearised acoustic theory, share a Maxwell-like relation [46,47].The time-harmonic equations are, where the constants β and ρ are the acoustic medium's compressibility and mass density respectively.Compared to photons, acoustic phonons, which are spin-0 quanta in this regime, do not give acoustic fields as rich a vector structure as light.Constraining the velocity vector is the longidutinality condition, that is ∇ × v = 0, a more restrictive analogy to light's transversality condition due to Gauss' law.Yet, v can still rotate and generate acoustic SAM, which is expressed (time-averaged) by, In time-harmonic acoustic fields, where the velocity v and displacement field r vectors are related by v = −iωr, particles in the acoustic medium have an elliptical motion.Meanwhile, an acoustic analogy to the Poynting vector can be defined by mixing the acoustic pressure and velocity fields, Taking the curl of P ac and following with use of Eq. ( 31), the acoustic SAM emerges in what appears to be the acoustic analogy to the presented electromagnetic spin decomposition: Due to the longitudinality condition ∇ × v = 0, the first term in Eq. ( 34) vanishes leaving the acoustic spin entirely proportional to the curl of the acoustic Poynting vector.Transverse phonons for which ∇ × v ̸ = 0 can occur in viscous fluids or solids, although in these media, the acoustic field could no longer be described by a linearised theory.Compared to the electromagnetic spin decompostion, however, Eq. ( 34) highlights the structural distinction between spin-0 and spin-1 fields as only vorticity in the flow of energy can generate SAM in a linear acoustic field.

On units
Discussions in this section are related to the 4-vector extension of the decomposition presented in section 4.2 of the manuscript.The field strength tensor F µν and conjugate field strength tensor G µν are not of the same units.The former has SI units kg s −4 A −1 whilst the latter has SI units m −1 kg s −5 A −1 .This comes from the form of the conjugate field strength tensor that can be derived using the companion Mathematica notebook [42]: Finally, we summarize all the SI units of the mentioned tensor and constants: Table 3: Summary of the SI units of every tensors and constants mentioned in this work

On the choice of gauge
The electric field and magnetic field, respectively E and H, can be expressed in terms of a scalar potential ϕ and a vector potential A, 36)

Figure 1 :
Figure1: The component parts of light's total spin angular momentum density, visualised at a single point (red circle) in the 3D interference of monochromatic plane waves.Enlarged electric and magnetic field ellipses and coloured vector arrows, representing the total spin S, electric and magnetic spins S e and S m , and the canonical s c and Poynting spins s p of the decomposition Eq. (4), are plotted.The electric and magnetic spin vectors are normal to the E and H polarisation ellipses, according to the right-hand rule and the sense of rotation of the instantaneous field vectors (see separated diagrams below the main combined image).Projected on the xy plane are the Poynting vector streamlines, which in this case organise a clear curling structure.This provides a visual aid for Poynting spin (large black arrow), which is informed by three-dimensional curling of the Poynting vector.In structured light where the E and H ellipses are de-coupled, the canonical spin and therefore chiral pressure (Eq.(5)) can point in a completely different direction to the any of the total spin or its electric and magnetic contributions.

Table 2 :
Comparison between acoustics, electromagnetism and linearised gravity.Table inspired by